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Abstract 

We study the full spectrum of spherically symmetric solutions in the five dimen- 
sional non-projectable Horava-Lifshitz type gravity theories. For appropriate ranges 
of the coupling parameters, we have found several classes of solutions which are 
characterized by an AdS^, dS^ or flat large distance asymptotic behaviour, plus the 
standard tail of the usual five-dimensional Schwarzschild black holes. In addition 
we have found solutions with an unconventional short or large distance behaviour, 
and for a special range of the coupling parameters solutions which coincide with black 
hole solutions of conventional relativistic five-dimensional Gauss-Bonnet gravity. 
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1 Introduction 

A novel quantum gravity model, which claims power-counting renormalizability, has 
been formulated recently by Horava [1]. This scenario is based on an anisotropy between 
space and time coordinates, which is expressed via the scalings t — )■ b^t and x — )■ bx, where 
2; is a dynamical critical exponent. For z 1 the UV behaviour of the model is governed by 
a non-standard Lifshitz fixed point, while for z = 1 we recover the well-known free Gaussian 
fixed point. In the Horava model, for three spatial dimensions, the suitable choice is z = 3. 

It is worth noting that in the Horava-Lifshitz (HL) gravity the four-dimensional diffeo- 
morphism invariance of general relativity is sacrificed in order to achieve power-counting 
renormalizability. The action of the model can be split into a kinetic plus a potential term, 
which both respect a restricted (3 -|- 1) diffeomorphism invariance. The interesting feature 
is that the kinetic term contains only second order time derivatives, while the potential 
term consists of higher order spatial derivatives of the metric components. This particular 
structure improves significantly the UV properties of the graviton propagator, and renders 
the model power-counting renormalizable. Moreover, in this way we avoid ghost modes 
which are usual in conventional higher order gravity models. 

For the construction of the potential term, there has been proposed llj the so called 
"detailed balance principle", which is inspired by condensed matter physics. The main 
advantage of this approach is the restriction of the large number of arbitrary couplings 
that appear in the bare action of the model. However, the physical motivation for a 
consideration such as the "detailed balance" is not clear [21 |3]. An alternative way for 
constructing an action is to include all possible operators which are compatible with the 
renormalizability of the model; this implies that all operators with dimension less or equal 
to six are allowed in the action (for the exact form of the action see [3]). 

As it has been already mentioned, the HL gravity violates local Lorentz invariance in the 
UV, however it is expected that general relativity is recovered in the IR limit. This implies 
a very special renormalization group flow for the couplings of the model. In particular, 
the parameter A in the kinetic term of the action (which measures the departure from the 
Lorentz invariance) should flow to unity, while the higher order couplings should vanish, or 
they should become appropriately small, in the IR. Note that, even though phenomenology 
suggests a particular IR limit there is no theoretical study supporting this behaviour. 

In addition, there are several other possible inconsistencies in HL gravity which have 
been discussed in several works. In particular, the absence of full diffeomorphism invari- 
ance introduces an additional scalar mode which can lead to strong coupling problems or 
instabilities (see for example |l]-[9] and references therein). However, these problems will 
not be addressed in the present paper. 

Apart from these problems, the HL gravity is an interesting quantum gravity theory, 
which has stimulated an extended research on cosmology and black hole solutions [TU] - 
[25] . In addition to general relativity studies, quantum field theory models in fiat space-time 
with anisotropy have also been considered [26j-[36j. 

Before proceeding it is important to mention that the HL gravity can be separated 
into two versions which are known as projectable and non-projectable. In the projectable 
version the lapse function N"^ depends only on the time coordinate, while in the non- 
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projectable version A^^ is a function of both the space and time coordinates. Although in 
General Relativity the projectable and non-projectable ansatz for the metric are equivalent, 
since they are connected via a diffeomorphism transformation, in the HL gravity the full 
diffeomorphism invariance is broken and they lead to two distinct theories. 

Issues connected with broken Lorentz invariance were studied in the spherically symmet- 
ric solutions of 4D HL gravity. In the case of detailed balance, such spherically symmetric 
solutions were found [11], but they exhibited an unconventional large distance asymptotic 
behaviour. The correct Schwarzchild-fiat asymptotic behaviour can be recovered if the 
detailed balance action is modified in the IR by a term proportional to Ricci scalar, and 
the cosmological constant term is considered to be zero ^2j. A similar study, in the case 
of non- vanishing cosmological constant, has also been carried out [13]. A generalization to 
topological black holes was obtained in [H]. Finally, a systematic study of static spheri- 
cally symmetric solutions of 4D HL gravity was presented in [15j where the most general 
spherically symmetric solution for A 7^ 1 and general coupling parameters was obtained. 

In this work we present a full study of spherically symmetric solutions in the non- 
projectable version of the five- dimensional Horava-Lifshitz gravity, for z = 4. For the 
construction of the 5D action we do not use the "detailed balance principle", but we 
include all the terms which are compatible with the renormalizability of the model. In 
particular, we can include all spatial curvature terms with dimension less than or equal to 
eight. However, the large number of possible terms, which are allowed in the action, leads 
to equation of motion of great complexity. For this reason we restrict our study only to 
terms of up to second order in the curvature. Also, we suppose that in the IR limit 5D the 
HL gravity reduces to the 5D General Relativity plus a bulk cosmological constant. A class 
of spherically symmetric solutions of the 5D HL gravity has been considered previously 
[17] . but only for a very specific choice of the couplings. 

Our main motivation in considering static solutions in 5D HL gravity is to investigate 
whether the rich spectrum of black hole solutions found in 4D (see Ref. [15]) also persists 
in 5D. It seems that the known static solutions of the HL gravity in 4D with A 7^ 1 do 
not have any obvious relation with the corresponding static solutions of the relativistic 4D 
gravity. In 5D however we found that there is a class of spherically symmetric solutions 
which after a proper identification of coupling parameters coincide with the known black 
hole solutions of conventional relativistic 5D Gauss-Bonnet gravity. 

Static solutions of 5D Gauss-Bonnet theory are well known Among them there 

is a black hole solution which has two branches (for a review see [3S])- The first branch 
is referred to as the Einstein branch while the second as the Gauss-Bonnet branch. Both 
branches coincide in the Chern-Simons limit. As we will discuss in the following, we find 
both branches of solutions and in addition these solutions can also be obtained for a dif- 
ferent combination of coupling parameters of the quadratic curvature terms than the usual 
combination that appears in the relativistic Gauss-Bonnet theory. We also find the black 
hole solution corresponding to the Chern-Simons limit with a particular choice of coupling 
parameters. This solution has also been found in [17] using the "detailed balance principle" . 

The paper is organized as follows. In section 2 we write down the action of 5D HL 
gravity. In section 3 we derive the equations of motion. In section 4 we analyze the static 
spherically symmetric solutions for a special choice of coupling parameters. In section 5 we 



-4- 



study the most general static spherically symmetric solutions of 5D Horava-Lifshitz gravity 
and finally section 6 contains our conclusions. 

2 5D Horava-Lifshitz gravity models 

In this section we introduce the notation for the so-called Horava gravity models in the 
case of four spatial dimensions (d = 4). These models are characterized by an anisotropy 
between space and time dimensions 

[t] = -z, [x] = -1 , (2.1) 

where z is an integer dynamical exponent. In order to derive the action of the model, it is 
useful to express the space-time metric in the ADM form 

ds'^ = -c^N^df + Qij {dx' - N'dt) {dx^ - N^dt) , (2.2) 

where c is the velocity of light, with dimension [c] = 2; — 1, and spatial components dx^/dt 
{i = 1,2,3,4). In addition, and Ni are the "lapse" and "shift" functions which are used 
in general relativity in order to split the space-time dimensions, and Qij is the spatial metric 
of signature (+,+,+,+). For the dimensions of "lapse" and "shift" functions we obtain 

[N]=0, [N,] = z-1. (2.3) 

In this paper the dynamical exponent z is set equal to 4. The 5D action of the model is 
constructed from a kinetic plus a potential term according to the equation 

S = [ dtd'^x./lilN {Ck + Cv} (2.4) 

in which d {D = d+1 = 5) is the spatial dimension and is the five dimensional Newton 
constant. 

We stress that the main motivation for considering models of this type is the con- 
struction of a power-counting renormalizable gravity model. However, in order to achieve 
renormalizibility, and simultaneously keep the time derivatives up to second order, we have 
to sacrifice the standard 5D diffeomorphism invariance of general relativity, which is now 
restricted to the transformation 

x' = x{x^,t), i=i{t) . (2.5) 

The kinetic part in the above Lagrangian of Eq. (12. 4p can be expressed via the extrinsic 
curvature as: 

Ck = {K''K,,-XK^), K,, = ^{~dtg^j + ViN, + V,N,}, z,j = l,2,3,4 (2.6) 

which is invariant under the transformations of Eq. ( 12. 5p . For the construction of the 
potential term we will not follow the standard detailed balance principle, but we will use 
the more general approach [21 13] , according to which the potential term is constructed by 
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including all possible renormalizable operators □, that have dimension smaller or equal to 
eight, hence we write 

t^V = C-R + ^R'^ + ^R? + ^AR^ + ^Ri + ^AR3 + ^A^R^ ■ (2-7) 

where the symbol A is defined as A = didi {i = 1, 2, 3, 4). 

The dimensions of the various terms in the Lagrangian read 

[R] = 2, [R^] = 4, [R^] = [AR^] = 6, [R^] = [AR^] = [A^R^] = 8 . (2.8) 

In this work we are mainly interested in the lowest order operator and the operator 
C^2, which contains contributions of second order in the curvature: 

'Cr = VOa + VlaR, Cn2 = r]2aR^ + V2bR''Rij + V2cR'''' Rijkl , (2.9) 

where we have used the notation R, Rij and Rijki for the Ricci scalar, the Ricci and the 
Riemann tensors = 1, 2, 3, 4), which correspond to the spatial 4D metric Qij. Note that 
in the case of three spatial dimensions the term R^^^^Rijki is absent, as the Weyl tensor 
in three dimensions automatically vanishes. However, in four spatial dimensions this term 
cannot be omitted from the action. 

The first term Cr is necessary in order to recover 5D general relativity with a cosmo- 
logical constant in the IR limit. The second term Cr2, includes all possible quadratic terms 
in curvature, and becomes important in the short distance regime of the theory. Moreover, 
?7oa plays the role of the cosmological constant, while r]ia, rj2a, V2b, and ?72c are dimensionful 
coupling constants with dimensions 

M=6, [rj2a] = [mb] = [V2c]=^. (2.10) 

In the present analysis we ignore higher order Lagrangian terms, of dimension six and 
eight. Although, we have not derived the detailed expression for the Lagrangian, it is 
worth writing some of the higher order curvature terms here, 

£^3 =R^ + RijR^R + £a/?2 = R^R + R'^ARij + (2.11) 
£^4 = i?^ + {RijR'^f + C^R3 = R^AR +..., C/^2R2 = RA^R + ... (2.12) 

where A = V^V*. The short distance effects of these terms may be important, but this 
topic is left for a future investigation. 

If this model is to make sense, it is necessary that the 5D general relativity (with 
a cosmological constant in our case) is recovered in the IR limit. Although there is no 
theoretical proof for this difficult question, we will assume that the renormalization group 
fiow towards the IR leads the parameter A to the value one (A = 1), hence 5D general 
relativity is recovered. Also, to obtain the Einstein-Hilbert action 

Seh = j dx^d^x^N [k^^k^^ -k^ + R + r]oa) , (2.13) 

we have to set rjia = c^, and 

A-« = ^{-9o«, + V,(f)+V,(^)} . (2.14) 

where the time-like coordinate xq is defined as xq = ct. 
^We have ignored the terms which violate parity, see also [3]. 
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3 Equations of motion 



We are looking for 5D spherically symmetric solutions of the Horava-type gravity model 
we constructed in the previous section. We use the following ansatz for the metric 



-N{rfde + f-\r) dr^ + r^d^l , 



(3.1) 



in which r is a radius coordinate that corresponds to the extra dimension, and dQl is 
the metric of a 3D maximally symmetric space, where k is the spatial curvature of 3D 
hypersurfaces and for k = 1,-1,0 we have a sphere, hyperboloid or 3D torus topology 
correspondingly. In what follows it is convenient to perform the transformation 



/(r) = k + r^Z{r) . 
Then the action of the model to second order in curvature terms is 



(3.2) 



iOTTCrs J 

(3.3) 

which can be put into the form 



S 



iV(r),Z(r), 



dZ{r) 
dr 



+00 



dr L 







N{r),Z{r) 



dZ{r) 
dr 



(3.4) 



after we integrate out the angular coordinates, where 

2 



dZ 
dr 



3'^ p(.f).M(.)(4).g,z) 



(3.5) 



and the coefficients P, M{Z) and Q{Z) are defined by the equations 

P = 3(3?72a + V2b + V2c) , 

M{Z) = 6{12r]2a + Sr]2b + 2r]2c)Z - 3r]ia , 

Q{Z) = 12{12r]2a + S7]2b + 2r]2c)Z^ - 12r/i,Z + r]oa • 



If we set 



V = i^V2a + V2b + V2c) , Q = (12r72a + ^'n2b + 2r72c) 



(3.6) 



(3.7) 



we can reduce the number of the free parameters of the model. This is possible because 
of the spherical ansatz for the metric, as it is given by Eq. (13.11) . In we set rjia = 1, by 
choosing a coordinate system in which c = 1, we obtain the simplified expressions 



P = 3r], 

M{Z) = 6gZ ~3 , 

QiZ) = 12qZ^ - 12Z + 7]aa . 



(3.8) 



^There is a more general ansatz for the metric, of the form ds^ — ~N{r)^dt^ + f ^(r) (dr + N"^ {r)dt)^ + 
r^dT,^ with nonzero shift N'^{r), but we have set N^{r) = [15] to simphfy the equations. 
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We note that the coefficient P is independent from r, while the functions M{Z) and Q{Z) 
do not depend exphcitly on the radius r. The Euler-Lagrange equations for the action (13.41) 
are 

d r dL\ dL _ d r dL\ _dL _ 
d^\dN')~dN~ ' d^\dZ')^dZ~ ■ 

The ffist equation of motion, the one for the function Z{r), reads 

p{r§)\M(Z)[r§)^QiZ)=0. (3.10) 

If we algebraically solve the above equation we obtain the ffist order differential equations 

(3.11) 

where H{Z) are the solutions of second order algebraic equation (l3.1Qp . For P 7^ : 

(3.12) 



rf^mz) 



H{Z) 



-M{Z) + (x^M{Zy -APQ{Z) 



2P 



where a is a sign; For P = we have only one solution: 



H{Z) 



Q{z) 

M{Z) 



(3.13) 



Now we can derive the equation of motion for the function N[r). If we set 



N{r) 



lN{r) 



fir) ' 

we obtain from the second Euler-Lagrange equation in (13. 9p : 



(3.14) 



dN{r) 
dr 



+ C{r)N{r) = 0, C{r) 



1 d{r^Gi) 



r^Gi dr 



rGi 



where 



G, = 2P[ rf] + M{Z) 



G2 = M'{Z) (^r^^+Q\Z). 
Changing variables from r to Z, Eq. (I3.15P becomes: 



(3.15) 

(3.16) 
(3.17) 



dN{Z) 
dZ 



C{Z)N{Z) = , C{Z) 



H{Z) 



4-a 



1 dGi 
Wi'dZ 



(3.18) 
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where 

Gi{Z) = 2PH{Z) + M{Z) , 

G2{Z) = M'{Z)H{Z) + Q'{Z) . (3.19) 

Finally, we emphasize that the parameter A does not appear in the equations of motion, 
so they only depend on three parameters: t], g and the cosmological constant ?7oa- This is 
similar to the 4D case [TB]. The reason is that we are looking for static solutions, hence 
extrinsic curvature terms do not contribute to the equations of motion, as they contain 
only time derivatives of the metric components. In addition, note that the parameter A 
appears only in the extrinsic curvature part of the action. Therefore, the solutions we will 
obtain in the following sections will be valid for arbitrary values of A. 



4 Spherically symmetric solutions, special cases 

In this section we study static spherically symmetric solutions for three special cases of 
the free coupling parameters t] and g: a.) t] = and ^ = 0, b) ?7 = and ^ 7^ 0, c) ^ = 
and rj 0. 



4.1 No quadratic terms, 77 = and g = 

If we set rj = g = in Eq. fl3.10p we find 

rf = f - 4Z . (4,1) 

from which it follows that 

-3Z+^ + % = 0, (4.2) 
or equivalently we take the simple solution 



f{r) = k + r'Z = k + ^r' + ^ , (4.3) 
J V ; r 12 3r2 ' ^ ^ 



where is a constant of integration. If we set 



C 

Kff = -Voa , = ' (4-4) 



the above equation takes the well-known form 



m = k - - ^ , (4.5) 



which is the standard AdS^ (for A^fj < 0) or dS^ (for Ag// > 0) or asymptotically flat (for 
Ag// = 0) Schwarzschild black hole solution of 5D general relativity with a cosmological 
constant. 
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4.2 r] = and g 

If 1] = Eq. fl3.8p implies P = 0, hence Eq. f l3.10p can be written as 

dZ _ Q{Z) 
~ ~ M{Z) ' 

In this case the function Q{Z) and M{Z) can be written as 

M{Z) = -3 + 6gZ , 
Q{Z) = r]oa - 12Z + 12gZ^ , 

and Eq. (14. 6p becomes: 

dZ _ 7]oa - VIZ + 12^Z2 
dr 

Integration of this equation yields: 



-3 + 



4 



(4.6) 

(4.7) 
(4.8) 

(4.9) 
(4.10) 



where C*^ is an integration constant which is related to the mass of the black hole. The 
algebraic equation ( I4.10p gives two solutions 



1 v/3(3 - ^r7o„)r4 - 12f?C^ 
Zir) = h a- , 

where cr is a sign {a = ±1), so for the function /(r) = k + r'^Z we obtain 



(4.11) 



f{r) = k + 



2g 



(4.12) 



In what follows we will assume that gCfj, < 0, because for gCfj, > the range of radius r 
has a lower bound (r > rmin)- This case will be discussed further in section 5.2. 

From the above equation fl4.12p we can extract the large distance asymptotic behaviour 
for /(r), which reads 



f{r) 



a 



+ 



1 



(4.13) 



Note that /(r) has a large distance limit only if 3 > grioa, or else there is an upper bound 
for the radius r. The asymptotic formula fl4.13p is of the form 



fir) 



12 r2 



(4.14) 



with 
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Depending on the values of the free parameters q and ?7oa, the asymptotic behavior is either 
AdS^ (for Keff < 0), or dS^ (for Ag// > 0), or flat (for Ag// = 0). Also, note that [/x] = 2. 
and Ae// is an effective 5D cosmological constant. The mass parameter of the black hole, 
when 7] = ^) = 0, is m = (87rG5)~^A~^^j/i. 

The Euler-Lagrange equations for N[r) yield 



+ - 



1 dGi 



Gi dZ 



(4.16) 



where 



M(Z) 



M(Z) 



Q\z) 



Then we obtain 



where we have taken into account equations (14. 7p . (14. 8 p for M{Z) and Q{Z). Finally, we 
find 

NirY = fir) . (4.17) 



4.3 Comparing with the 5D Gauss-Bonnet gravity 



It is worth noting that the spherically symmetric solutions we obtained in the previous 
section, for rj = and g ^ 0, are identical with the corresponding solutions of the 5D 
relativistic Gauss-Bonnet (GB) gravity. The action of the GB gravity is given by: 



S = ^Jd^x^\{R^^^-2A + dG}, D = 5, 
where the GB density G is 

Q ^ j^i5)abcdj^(^)^^ _ 4^Jl(5)abj^(J>) ^ (^^(5)^2 ^ a, b, C, d = 0, 1, . . A 



(4.18) 



(4.19) 



Note that in the GB gravity, the definition of the symbols R^^\ and i?(5)«^'cd based 

on the relativistic 5D metric g^^J , where a,b,c,d = 0,1,. ..4. In addition, d is the Gauss 
Bonnet coupling and A is the 5D cosmological constant. The static spherically symmetric 
solutions of the GB gravity in AdS space, for D = 5, are of the form 



ds^ = -f{r)de + r\r) dr^ + r^dQl 



and 



fir) 



4a 



a 



±1 



(4.20) 



(4.21) 



in which fig is a constant of integration which is related with the mass of the black hole, 
and the parameter = — 2A corresponds to a negative bulk cosmological constant. 
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If we replace 



a — )■ 



2 , VOa 



— 



a, 



(4.22) 



2 12 3 

in equation f l4.2ip . we recover the black hole solution of Eq. f l4.12p . for the specific case 
?7 = and 7^ of the previous section. 

Note, that the condition 77 = 3ri2a + V2b + V2c = is satisfied in the case of GB coefficient 
r]2a = 0,, rj2h = —4a and 7y2c = a, but there are other different combinations of the coupling 
parameters 7720? V2hi V2c which give rj = and g 0. This is a very interesting result which 
merits further investigation. Note also that the relation 1 — = corresponds to the 
Chern-Simons limit of GB gravity. 



4.4 ^ = and 77 7^ 
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Figure 1: A typical plot of /(r) and N{r)'^ versus r, in the case of g = and ?7 7^ 0, for k = 0, 
r]oa = — 1 (positive cosmological constant), = 10, r] = —0.1 (top) and r/ = +0.1 (bottom) 

If we assume that g = 0, Eqs. ( 13. Sp yield 

P = 3r], M{Z) = -3, Q{Z) = -12Z + 7]oa . (4.23) 
By solving Eq. ( 13.10p we obtain the following two solutions 



^dZ ^ 3 - gy/Q - 12r]r]oa + UAr]Z 
dr 6r] 
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with cr = ±1. The above differential equation can be integrated analytically, to give 

where is an integration constant. 

In this section we are interested only in solutions with cr = 1 and > 0. The other 
cases lack a short distance limit (for details see Appendix A) and hence their interpretation 
is problematic, as we explain in section 5.2 below. 

From Eq. f l4.25p we obtain 

'C„ 



- 12r]r]oa + U^VZ - I = Wl j , (4.26) 

Note that Wl{x) is the Lambert function, which is defined as the real solution of the 
equation e^^^^^Wiix) = x. We recall some of the properties of this equation: (a) for 
X < — 1/e it has no real solutions, b) for — 1/e < x < it has two real solutions and, 
c) for a; > it has a unique real solution. In the case b) we define the function Wl{x) 
by demanding that — 1 < Wl{x) < (the other set of solutions, which lies in the range 
(— oo, — 1), is not considered). 

Now, from Eq. (14.261) we find that 



so we find for the function f(r) = k + r Z(r) : 

/(r) = k + ^-^r' + ^ (^Wl + 2W, (^^]\. (4.28) 

The large r asymptotic behaviour of Eq. f l4.28p is found to be 

12 8T]r^ \f 

which is of the standard form 



f{r) = k+'-^r' + ^^+Oi-^) , (4.29) 



m^k-^r^-^^, Kj, = -Voa, /^ = -^- (4-30) 



Now, Eqs. dMD, f l338|) and i KWf yield 



where the function N{Z) satisfies the equation 

dN{Z) 



+CiZ)NiZ) = 0. (4.32) 
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Solving this equation and choosing the constant of integration appropriately, we get 

o [1^9-12^^+144772-1] 

NiZ) = , . (4.33) 

Note that for a = 1 and > 0, if we take into account Eqs. f l4.26p and f l4.33p . the function 
N{rY can be expressed in the following closed form 

N{rf = f{r)N{Z{r)f = — . . (4.34) 

In the large r regime we find, from the above equation, that: 

N{rf = f{r) (^1 + ^ + (^^)) , (4.35) 

hence in the large distance limit we recover the standard asymptotic behavior A^(r)^ ~ /('")• 
In Fig. [T]we give a typical plot of the functions /(r) and A^(r)^ for zero spatial curvature 
A; = and positive effective cosmological constant Ag// = 1. We see that /(r) is finite (in 
particular /(O) = 0) but N{rY blows up at the singularity = 0. We also see that 
only when the coupling rj is positive there exists a horizon, while for negative 77 there is a 
naked singularity. In the case of non-zero spatial curvature the horizon r/j is determined by 
solving equation f{rh) = —k, {k = ±1). It is possible to see by inspection of Fig. [H that 
for k = ibl, we have the same situation which holds for zero spatial curvature. Finally, we 
observe that /(r) and A^(r)^ tend rapidly to their common asymptotic behavior {dS^), as 
it is expected from the analysis above (see Eqs. f l4.29p and f l4.35p ). 

5 Static solutions in the generic case {rj ^ and Q ^ 0) 

In the generic case, 1] and g ^ 0, we obtain from Eqs. (13. 8p . fl3.1ip and f l3.12p 



dZ 3-6gZ + a^9-12r]oaV + ^Qi0-'iv)iQZ^-Z) 

r— = . (5.1) 

dr 6?7 



If we make the replacement 



in Eq. (15. ip . we find 

= H{y) , H{y) = -- (^gy + a^A + By^) , (5.3) 



ri 

where the new parameters A and B are defined through 



A = -3r]oaV + — , B = g{g-A7]) . (5.4) 
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For the computation of the function N{r), Eqs. 03. Sp . 03.18P and f l3.19p yield 

dN{y) 



dy 



where C{y) is given by the equation 



C{y) = - 



C{y)N{y) = 



B gy + ay/ A + By"^ 



(5.5) 







A + By^ 



(5.6) 



In the following sections we analyze two cases for the parameter B: 1) B > 0, 2) B < 0. 
The mathematical details for the derivation of the final formulae for /(r) and N[r) are 
given in Appendix B. 



5.1 B>0 

By integrating Eq. fl5.3p for i? > we obtain 



0y 



ay/A + By^ VBy + a^/ A + By^' 



(5.7) 



Note, that in this case the parameter A can be positive, negative or zero. For the function 
iV, we obtain from Eqs. (15. 6 P and (15. 5p : 



c 



N 



^A + By"^ y/By + a^/A + By^ 



g 



where and are constants of integration. Note, that must be always positive 
(C^ > 0). Also, the constant Cn is fixed if we demand — )■ 1 for r — )■ +oo, as we will 
discuss later in this section. An alternative expression for A^ can be found, if we take into 
account Eqs. (15. 7p and (15. 8p . In particular, we obtain that 



N{y) 



AT 



\gy + cta/A + By"^ 



C 



N 



CmC, 



(5.9) 



For B > 0, two cases for the ratio ^/B / g will be discussecE]: a) \y/B/g\ < 1 and b) 
\\/B/g\ > 1. The case B = will also be examined separately in section 5.3. 

5.1.1 \VB/g\ < 1 and A>0 



For \\/B/g\ < 1 and A > 0, we can verify that Eq. (I5.7p above, has two solutions for y {yi 
and 1/2) for a given value of the radius r in the range [0, +00). Hence, the function /(r) has 
two branches fi{r) and f2{f'), which are exhibited in the left part of Fig. [21 when a = 1 



•^Note that if \VB/g\ = 1 we have no solution at all. 
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for two typical values of ^ = ±10). However, only one of them has a horizon for zero 
spatial curvature {k = 0), while the other represents a spherically symmetric solution with 
a naked singularity. In addition, as we see in Fig. [21 for ^ > the black hole solution is 
AdS^ asymptotically, while for < is dSr, (note that g should be non-zero in the case we 
examine here). 

In what follows, we will assume that A > 0, as for negative A the solutions have no 
large distance limit when \y/B / g\ < 1, hence they will not be examined here. 
In particular for / g\ < 1, if we take into account Eq. (15. 7p . we obtain 



• r — 7- +00 ^ y — 7- yo (large distance asymptotic behaviour) 

• r — 7- ^ ?/ — 7- ±oo (short distance asymptotic behaviour) 

where the plus and minus signs above correspond to the two branches of the function f(r). 
In addition, 

yo = -sgn y^jAl (5.10) 
is the unique solution of the following equation 



gy, + a^A + Byl = Q . (5.11) 

If we expand Eq. (15. 7p around we can find the large distance asymptotic behavior for 
y{r) which reads 

y{r) . .0 ± ^ + O = |,| l^o(^^-gl^ c, , (5.12) 



yry'"^ 3{g^-B) 
hence 

/(r) = A: + r^Z = fc + r^(l-|)^/(r)^A:+(l-|)r^±ii + 0(l) , (5.13) 

which has the standard asymptotic behaviour, AdS^, dS^ or flat, depending on the values 
of the free parameters of the model. The plus and minus signs in the above asymptotic 
formulas give rise to the two branches of the function /(r). 

Now, for the function N{y), if take into account Eqs. (15. 9p and (I5.12p . we obtain the 
following large distance asymptotic behaviour 

iV(y(r))^l + 0(^l^ , (5.14) 

where the constant of integration Cn has been set to 

Cn = -(Tg yo\yo{VB - g)\ e , (5.15) 



■^For A < and \VB/g\ < 1 the left hand side of Eq. (|5.7|) is never zero, so the radius r has an upper 
bound. 
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in order to satisfy the condition N{yo) = 1. 

The asymptotic behaviour of Eq. fl5.14p is verified graphically in Fig. [31 We observe 
that the function A^(r)^ tends rapidly to its asymptotic behaviour which is identical with 
that of /(r), as it is given by Eq. f l5.13p . Also, we would like to note, that the above 
analysis is valid only when A > 0. For negative A we see that the left hand side of Eq. 
(15. 7p cannot vanish (see Eqs (I5.10p and (I5.1ip ), so the radius r has an upper bound. We 
conclude that this class of solutions lacks physical interest, since the function /(r) has no 
large distance limit, so it will not be examined here. The case A = is also examined 
separately in section 5.1.4. 




Figure 2: A typical plot of /(r) versus r, in the case oi k = 0, A = I, B = l,Cf^ = l,a = 1. In 
the left part of the figure the parameter g equals —10.0 (top), or +10.0 (bottom), while for the 
right part of the figure g = +0.5 (top) and g = —0.5 (bottom) 

5.1.2 > 1 and A > 

In the right part of Fig. [2] we have plotted the function /(r) when \\/B / g\ > 1 (cr = 1, 
Q = ±0.5). As we see, in contrast with the case \\/B/g\ < 1, now there is only one branch 
for the function /(r). However, there is an even more significant difference, as this class 
of solutions does not exhibit the standard AdS^, dS^ or fiat asymptotic behaviour, as it 
is shown in the following analysis. For negative a {a = —1) the behaviour of /(r) does 
not change significantly, so we do not give a figure for reasons of space. The asymptotic 
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Figure 3: A typical plot of the two branches of N(r)'^ versus r, in the case of > and 
\VB/q\ < 1, for k=0, A=l, B=l, = 1, a = 1 and q=10. We observe that N{ry ~ /(r) when 
r — +00, as expected. 

behaviour of /(r), for negative and positive a, is described in Eq. fl5.1.2p below, and /(r) 
obeys the same power law for both values of a. 

Although we will assume that A > 0, there are also solutions for negative A in the 
case \\/B/g\ > 1 . These solutions do not have any new features, as we comment in the 
following section 5.1.3. 

For \\/B / g\ > 1 and A> 0, from Eq. (15. 7p above we obtain 

• r — )• +00 ^ ?/ — 7- e(+oo) (large distance asymptotic behaviour) 

• r — 7- ^ ?/ — i- — e(+oo) (short distance asymptotic behaviour) 

where e is a sign defined as e = sgni^ga). Note, that the left hand side of Eq. (15.71) does 
not vanish for a finite value of y {y = yo), the way it did in the previous section. 

For the large distance asymptotic behaviour of y{r), if we take into account Eq. (15. 7p . 
we find: 

4|el 



?/ ^ Ce rv^-iei , (5.16) 

where the form of the coefficient depends on the sign of the parameter e = sgn{ga), 
according to the equations: 

(J = ( 1 A ^/s-|e| / VB+\g\ \ Vs-lel jf ^ ^ ]^ 
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The leading terms for /(r) are 



/(r)~A; + - '—r^^\ , (5.17) 

which is not the standard asymptotic behaviour for a 5D black hole solution, since it is 
proportional to r^+'', 5 = > ' rather than r^. The large distance asymptotic 

behaviour for the modified lapse function N{r), if we take into account Eqs. (15. 9 p and 
(EH, is 

Cat _4._i^M_ 

N{r) ~ ^ ^ r ^-1^1 . (5.18) 

Cf^y/ B\g + sgn(g)\' B\ 

We see that A^(r) — t- for large r, which implies a violation of the common large distance 
asymptotic behaviour A^(r)^ ~ /(^)- Finally, from the above equations (I5.17P and (I5.18P 
we can determine the leading term for the lapse function N{r), which reads 

N{r) ~ — — r ^-lel . (5.19) 

3vi?|f) + sgnlg)^ B\ 

We observe that also the lapse function N{r) vanishes for large values of r. 



5.1.3 \VB/g\ > 1 and A<0 

In fact, the solutions for \\/B / g\ > 1 and A < 0, is a mixture of solutions that were 
presented is the previous two sections. More specifically every solution consists of two 
branches. One of them has exactly the same short and large distance asymptotic behaviour 
with the solutions of section 5.1.1. This is mainly due to the fact that the left hand side 
of Eq. (15. 7p vanishes for a finite value of y. The other branch is similar to the solutions 
of section 5.1.2 but it lacks a short distance limit. We will not discuss these cases further 
since their behaviour is similar to the solutions we have already discussed. 



5.1.4 B >0 and A = 

If we set A = (or r]oaQ = 3) in Eq. (15. 7p . we find 



4e 



y[r) = C,ri^-^), Co = a{-j=] ' (5-20) 



2v^ ) \ C, 



so we obtain for the functions /(r) and N{r) 



f{r) = k + --^r'^~^ (5.21) 



^We remind the reader, that we can set A — only when B is positive or zero. 
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and 

Cat " v^+e 

N{r) = /75i ^ • (5.22) 

C^VB\g + VB\ 

As we see the above solutions exhibit an unconventional asymptotic behaviour which is not 
of the type AdS^, dS^ or flat 5D Schwarzschild form. 



5.2 B <0 




Figure 4: /(r) versus r, in the case oi B < 0, for k=0, A=l, B=-l, (7^ = 1, cr = 1 and q=-5 
(top), +5 (bottom) 



For i? < 0, if we integrate Eq. (15. 3p . we find 



+ a^/A-\B\ 



— - tan ^ ' — ^ ' 



^4 



(5.23) 



Note that in this case the parameter A should be always positive {A > 0). Eqs. (15.51) and 
(ESD yield: 



N{y) 



N 



'\B\y 



^A-\B\y^ 



(5.24) 



^A-\B\y^ 

Now, if we use Eqs. (I5.23p and (I5.24p we can get an alternative expression for according 
to the equation 

Cat 



(5.25) 
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Figure 5: N{r)'^/ f{r) versus r, in the case of B < 0, for k=0, A=l, B=-l, = 1, a = 1 and 
£•=+5. 

If we take into account the restriction from the square root in Eq. fl5.23l) . we find that y 
varies in the finite interval \y\ < i/|4|. 



In Fig. Hlwe have plotted the function /(r), for positive a (a = 1), two typical values 
of g {g = 5 > and g = —5 < 0) and zero spatial curvature k = 0. As we see /(r) consists 
of two distinct branches /i(r) and f2{r), with a dS^ (for ^ < 0) and AdS^ (for ^ > 0) 
large distance asymptotic behaviours correspondingly. In contrast with the previous case 
{B > 0) the range of radius r terminates at a lower non-zero bound r^mi (for /i(^)) and 
fmin2 (for f^i'f'))- FoT zcro Spatial curvature, A; = 0, we observe that for g < both the 
two branches have a horizon, while for g > we do not obtain horizons at all. Notice that 
the lower branch in the bottom of Fig. Hlmay have a horizon if A; = —1. In Fig. Owe see 
that the function N{rY/ f{r) tends to unity as r tends to infinity, but N{r)'^/ f{r) blows 
up when r — )■ rmini,2- We observe that in this case there is an infinite set of singularity 
points which lie on a 4D hypersurface of constant radius rmini,2- (This situation is to be 
compared with the standard case of a black hole solution where the singularity point is 
located at the origin). These solutions may have a physical relevance in the case where the 
singular shell is protected by a horizon, for example the branches of /(r) with negative g 
develop a horizon as we pointed out above. 

Also, in Fig. |6l we have plotted the function /(r) for negative a [a = —1) for the 
same two typical values of g and A; = . Note that there are no significant differences if 
we compare with the corresponding figure for cr = 1 (Fig. Hj). However, we see that for 
negative a only one of the two branches of /(r) possesses a horizon, while the other exhibits 
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Figure 6: /(r) versus r, in the case of B < 0, for k=0, A=l, B=-l, (5^ = 1, cr = — 1 and ^>=-5.5 

a naked singularity (in fact we have a naked singular shell). 

Although in the above mentioned figures we have used specific values for the free pa- 
rameters, our conclusions are quite general as the values we chose represent a wide range 
of the parameter space. 

In what follows we try to better understand some of the properties of the solutions by 
using the analytical formulae of Eqs. fl5.23p and fl5.24p above. First, we would like to stress 
that the transcendental equation fl5.23p gives two solutions for y for a fixed value of the 
radius r, as it may also be seen in Figs. H] and [51 In particular the two branches fi{r) 
and f2ir) correspond to the two intervals of the parameter y: 1) [— |^,2/o] and 2) [yo, -^]. 
Now, if we take into account Eq. fl5.23p we find that for large r the function y{r) tends to 
a constant value, or equivalently that 

r — )■ +00 =^ y ^ yo . (5.26) 

We note that 

is the unique solution of the equation 

gyo + a^JA-\B\yi = . (5.28) 

Note, that for negative B {B < 0) there is no short distance limit with r — )■ 0. We find 
that r > rmin, where r^m is a non-zero lower bound for the radius r, which is given by the 
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equation 

-1/4 



e^T. (5.29) 



It turns out that: 



lim t/ = ±^/A . (5.30) 



The large distance asymptotic behaviour for y{r) is obtained by expanding Eq. fl5.23p 
around yo : 



B\ . _W VIS 
tan 



l/W-?/oT^ + 0(l), /i= ,,'f'^ro. ^ ' (5-31) 



r4 ' V^V 3(^2 + 

hence for the function /(r) we find the corresponding asymptotic behaviour 

which has the standard AdS^, dS^ or fiat asymptotic behaviour, depending on the values 
of the free parameters of the model. 

For the function N{y{r)), when the radius tends to infinity, we find 

N{y{r))c^l + o(^-^^ , (5.33) 

where the constant of integration Cat has been set to 

Cn = -(ygyoe f ^ ^ , (5.34) 
to satisfy the condition N{yo) = 1. 

5.3 B = 

Although the special case of i? = is included in section 5.1.1, we present the relevant 
results separately here, since /(r) can be expressed explicitly as a function of r. If 5 = 
and A > 0, from Eq. ( 15. ip we obtain 

= -I (gy + aVA) . (5.35) 



dr 7] 

From B = g{g — 4i]) = and g 0, we conclude that g = Arj, hence the above equation 
f l5.35p can be written as 

= -- (^gy + aVTj , (5.36) 



Q 

which yields 



T ^ , (5.37) 



g gr"^ 



23 



while for the function /(r) we find 



3^ 



) 




(5.38) 



This has the standard form of a AdS^, dS^ or fiat solution. The above formulae can also 
be obtained if we set i? = in Eqs. f lS.lOp . f l5.12p and f l5.13p of section 5.1.1. Finally, from 
Eqs. (15. 5p and (15. 6p . choosing suitably the integration constant, we find that 




(5.39) 



It is interesting to note that (at least the positive branch of) the solution (I5.38P coincides, 
after a proper identification of the parameters, with the black hole solutions found in 
Lanczos-Lovelock theories [39]. The Lanczos- Lovelock action is a polynomial of degree [D/2] 
(the integer part of D/2) in the curvature and in 5D it has solutions of the Schwarzschild- 
AdS type. 

6 Conclusions and discussion 

We studied static spherically symmetric solutions in the framework of the 5D Horava- 
Lifshitz gravity. We considered an action consisting of terms of up to second order in the 
curvature and we solve the theory with a non-projectable spherically symmetric ansatz for 
the metric. The black hole spectrum we found is controlled by three parameters t], g and 
riQa, where ?7oa is a cosmological constant. The black hole solutions we found do not depend 
on the parameter A which measures the departure from Lorentz invariance as it appears 
only in the extrinsic curvature part of the action. 

We presented a full analysis of 5D Horava-Lifshitz static solutions scanning the values 
of the free parameters i], g and rjoa, which can be positive, negative or zero. This analysis 
comes as a generalization and extension of the 4D case studied in [T5j. More specifically, 
we obtained three main sets of solutions: the two special cases (r/ = 0, g ^ 0), and (r/ 7^ 
0, ^ = 0) and the generic case [r] ^ 0, g ^ 0). In all cases we obtained analytic black hole 
solutions which have the standard AdS^, dS^ of fiat asymptotic behaviour, plus the well- 
known 1/r^ tail. However, we also obtained solutions with an unconventional short and 
large distance asymptotic behaviour. For example, in the generic case, we found solutions 
with an asymptotic fall-off which is stronger than the AdS^ or dS^ asymptotic behaviour. 
We have also found solutions (in the cases with 77 7^ 0) in which the asymptotic behaviour 
is the usual one, but the radius has a lower bound r^m different from zero. Also, in many 
cases we obtained solutions with a naked singularity. 

We also found static solutions which, after a proper identification of coupling parame- 
ters, coincide with static black hole solutions of relativistic gravity theories with quadratic 
curvature correction terms. One class of these solutions consists of the Schwarzschild-AdS 
black hole solutions of five- dimensional Lanczos-Lovelock gravity theories. Another class 
of solutions contains the well-known Gauss-Bonnet black hole solutions. The interesting 
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result we obtained in our investigation is that the non-relativistic solutions of the HL grav- 
ity corresponding to the Gauss-Bonnet solutions can be obtained for various combinations 
of the coupling parameters t] and g and not just the standard Gauss-Bonnet combination. 
This may be attributed to the fact that the HL static solutions are insensitive to the cou- 
pling parameter A, so they hold even if A 7^ 1 (the value which signals the breaking of 
Lorentz invariance). 

We do not have a full understanding of the quantum UV structure of the Horava-Lifshitz 
gravity. The /3 functions for the UV marginal couplings have not yet been calculated, so 
the claim that in the IR the A = 1 is a fixed point is an assumption. The fact that a class 
of our solutions coincides with relativistic Gauss-Bonnet black hole solutions is suggesting 
that in the running of couplings towards the IR, the Gauss-Bonnet regime is reached before 
the 5D gravity is attained. After all the Gauss-Bonnet theory is an UV correction of 5D 
gravity. 

Important issues remain to be investigated. One of them is the stability of our solutions. 
For example, the stability of the class of solutions we found which coincides with the Gauss- 
Bonnet solutions is an interesting issue to be studied. The stabihty of the Gauss-Bonnet 
static solutions has been extensively studied |30]. It was found in [H] that one branch of 
these solutions suffers from ghost-like instability up to the strongly coupled Chern-Simons 
limit where linear perturbation theory breaks down. In our case because the Lorentz 
invariance is broken in the UV this behaviour could be different. Therefore, a careful 
stability analysis is required. 

Other interesting issues are the thermodynamic properties of our solutions, or the con- 
tribution of terms of higher order in the curvature, which have been omitted in the present 
work. It would also be interesting to generalize the solutions we presented here in the 
presence of electric charge. 

Finally, one field that our findings can be applied is the extra-dimensional gravity 
theories, in particular the brane world models. Note that in contrast to the standard 
vacuum of Randall- Sundrum |12], the five-dimensional AdS black hole vacuum does not 
preserve 4D Lorentz invariance on the brane, which may have interesting phenomenological 
implications SI] • For appropriate ranges of the coupling parameters, we have obtained 
solutions with an AdS^ large distance asymptotic behaviour, plus the standard tail 
of a usual 5D Schwarzschild black hole. These solutions may serve as backgrounds in 
the framework of brane worlds models, but now there is additional advantage: the starting 
point is a renormalizable theory such as the 5D HL gravity, in contrast with the 5D General 
Relativity, which is non- renormalizable and requires a UV complement. 
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8 Appendix A: Asymptotic behaviour analysis in the case g = 
and 77 7^ 

In this appendix we will examine the short and large distance asymptotic behaviour of 
the solution (14.251) in the case g = and r] ^ 0. In particular we will examine two cases I) 
(7 = 1, and II) a = —1 for the equation 

(^1^9 _ ur^r^,^ + UAr^Z - l) ^if V^^-^^^^+i^^^^-i] = Si . (g.l) 

As we will see, the only interesting case is that for a = 1 and > 0. The problem in 
other cases is that Z{r) is not defined on the whole interval [0, +00). 



8.1 Case I {a = I) 

Although, the case a = 1 and > has been examined in detail in section 4.4, we 
summarize our results here 

• r-^+oo^Z-^^ (and Z > ^) 

• r — 7- ^ Z — i- sgn{ri){+oo) 

For > 0, we see that Z{r) is well defined in the range [0, +00). 
For (7 = 1 and < we find 

• r^+oo^Z->^ (and Z < ^) 

• r — 7- ^ the limit does not exist, as there is a lower bound r > rmin 

For Cfj, < 0, the function Z(r) is defined in a range of the form [r^in, +00), where r^m > 0. 
In both cases the large distance asymptotic behaviour is given by Eq. (14.291) in section 4.4, 
which is of the standard form 



fir) ^ k - - 4 . (8.2) 

12 



In addition, for < 0, the Tmin can be determined if we take into account that the 
Lambert function has a lower bound Wl{x) > — 1/e, we then obtain 



_ -eC^ , (8.3) 
and 

Z(r„.„) ^ . (8.4) 

In order to derive the above equation we used that Wl{—1) = —1/e. 
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8.2 Case II {a = -1) 

For cr = — 1, from Eq. (18. ip we see that the parameter should be negative, hence we 
examine only the case (7^1 < 

• r — )■ +00 ^ Z — 7- sgn{ri){+oo) 

• r — )■ ^ the hmit does not exist, as there is a lower bound r > rmin 

We see that for a = —1, the function Z(r) is defined in a range of the form [rmin, +oo), 
where rmin and Z(rmm) are given by Eqs. f l8.3p and (18 .40 above. In addition, the large 
distance asymptotic behaviour in this case can be estimated from Eq. (18. ip if we keep only 
the exponential term 

Z(r)^-l-ln|^l , (8.5) 



hence 



16?7 \ r 



/(r)^A; + r2^1nf^l , (8.6) 



while for A^, from Eq. (I4.33p . we obtain 



NiZir)) ^ . (8.7) 



We conclude that for cr = — 1 (and g = 0, r] ) the large distance asymptotic behaviour 
is not of the standard form (18. 2p . so we will not study this class of solutions further. 



9 Appendix B: Technical details in the generic case [f] ^ and 

In this section we present same mathematical details for the derivation of formulae of 
Eqs. (15. 7p and (I5.23p . In particular we have to integrate Eq. (15. ip 



^dZ ^ 3-60Z + gy/Q - 12r7o,r/ + 36(g - A7]){gZ^ -Z) 
dr 6r] 

In order to perform the above integral we make the replacement Z = ^ — |, then we obtain 



dr 



H{y) , H{y) = (^gy + a^/JTW) , (9-2) 



where 

9n 

A = -3voaV + — , B = g{g-4r]) . (9.3) 
Q 

Note that 

2(^2 _ 5) 8' 
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so we have — B = Aqt] ^ 0. For 5 7^ 0, we obtain 



c + ln(r) 



dy 



-V / dy 



~ In I {g^ - B)y'^ - A\+a7] / dy 



gy-a^A + By^ 
(^2 _ B)y^ - A 

VaTW 



(9.4) 



(9.5) 



(^2 _ B)y^ - A ' 

where c is an integration constant. In order to calculate the second integral in Eq. (19. 4p 
separately, we write the corresponding integrand in the form 



B 



+ 



Ag^ 



1 



ig^~B)y^-A g^ - B \^A + By^ ^ A + ByHio" - B)y^ - A) ^ 
For S > we obtain 

y/B 



(^2 _ 5)y2 _ ^ 



g'^-B 
Q 



In 



while for i? < we get 



dy 



^A-\B\y^ 
{g^ + \B\)y^-A 



2{g^-B) 

1 

Q 



In 



By+^A + RiP 
gy-^ATW 



gy+^A + B^ 

\B\y 



|B|tan 



-1 



^A-\B\y^ 



In 



gy~^ATW 



gy+^ATW 



where we have taken into account that 
Ag" 1 



Ag^ 



,/aTW - B)y' - ^) + By')'/' (^1 

and the well known relations 



A+By^ 



d 



A 



dy\^/ATW {A + By^y/^' J l-x^ 2 



dx 



In 



1 — X 



1 + x 



X 



Va + W 



Moreover, we will use the identity 



iln 

2 



gy-^/XTW' 



gy+^ATW 



-\\n\{g' - B)y' - A\+\yi gy+ ^A + By^ 



(9.6) 



(9.7) 



(9.^ 



(9.9) 



(9.10) 



(9.11) 



For B > from Eqs. (19. 4p . ( 19. 7p and (19. lip , if we take into account that g^ — B = Agrj, we 
obtain 



c + ln(r) 



a-1 



- I- 



In \ {g' - B)y' - A\ 
gy+^/JTW , 



aVB 



In 



By+^/XTW 



(9.12) 
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or equivalently 



d + In(r) = — In 
4 



B 



4g 



In 



By + a^/A + By^^ 



(9.13) 



where c' is a new constant which is defined appropriately. Finally, we obtain 



+ o^A + By'^ \fBy + A + 5?/2 



4 ' 



(9.14) 



where = e'^ is the final integration constant which is related to the mass of the black 
hole. For i? < 0, in a similar way, if we use Eq. (19. 8p instead of (19. 7p . we obtain 



QV 



+ a^A- 



^4 



Finally, for the computation of N{r) we use the following equations 

dN{y) 



dy 



where Ciy) is given by the equation 



C{y) 



-C{y)N{y) = Q 



B oy + cta/A + By'^ 



Q A + By^ 
and we find in a similar way, as it is described above, that 

N{y) = ^ 



^A + By^ ^fBy + a^A + By' 



for i? > 0, and 



N{y) 



C 



N 



tan 



\B\y 



^A-\B\y^ 

for i? < 0, where Cfq is an integration constant. 



y/A-\B\y'^ 



(9.15) 



(9.16) 



(9.17) 



(9.18) 



(9.19) 
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